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Abstract: The large-order behaviour of QCD is dominated by renormalons. On the other 
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f^ I hand renormalons do not occur in conformal theories, such as the one describing the infrared 

fixed-point of QCD at small (3q (the Banks-Zaks limit). Since the fixed-point has a pertur- 
bative realization, all-order perturbative relations exist between the conformal coefficients, 
which are renormalon-free, and the standard perturbative coefficients, which contain renor- 
malons. Therefore, an explicit cancellation of renormalons should occur in these relations. 



D , The absence of renormalons in the conformal limit can thus be seen as a constraint on the 

'^ ' structure of the QCD perturbative expansion. We show that the conformal constraint is 

non-trivial: a generic model for the large-order behaviour violates it. We also analyse a 
^ ■ specific example, based on a renormalon-type integral over the two-loop running-coupling, 

H , where the required cancellation does occur. 
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In QCD, as in other quantum field theories in four dimensions, the running- 
coupling depends logarithmically on the scale. As a consequence, the perturbative 
expansion is characterized by a renormalon factorial increase, which emerges from 
integrating over the running-coupling in loop diagrams [0-0|. 

Renormalons are believed to dominate the large-order behaviour of the series. 
In many physically interesting cases, their dominance sets in rather early 0]. Then 
renormalon resummation, i.e. the summation of the specific contributions associated 
with the running-coupling to all orders, becomes important for precision calculations 
(see e.g. !-§). 

The idea that running-coupling effects can be resummed to all orders is based on 



an analogy with the Abelian theory |T0, |TT|, where there is a systematic skeleton ex- 



pansion. In the absence of the latter, the separation of running-coupling contributions 
from the conformal part of the expansion is not well defined. In practice, in most 
applications, the resummation is restricted to the level of a single dressed gluon, 
where the (Abelian) \arge-Nf calculation is sufficient to determine the renormalon 
contributions. The resummed result can then be written as 

RoiQ') = f<k')<P{kVQ')^ (1) 

= a{Q') + r?^(3oaiQr + (rf /?„' + ri^^^) a{Qy 
+ (rf /5o^ + ^rf A/5o + rf)/?,) aiQY + • ■ • , 

where fc^ is the virtuality of the exchanged gluon, a(fc^) = as(A;^)/7r is the QCD 
running-coupling satisfying the renormalization group equation, 

^^ = Piaik')) = -Poaik'f - PMk'f - ^2a{k'f -■■■ (2) 

and (f) is an observable dependent Feynman integrand for a single gluon exchange 
diagram, which is interpreted as the gluon momentum distribution function 0. 

Equation (P, should be viewed as the leading term in a hypothetical skeleton 
expansion fllj: i?(Q2) = Rq{Q'^) + i?i(Q2) + . . .. it has also proved useful §, |, |] to 
view it as a model for the all-orders structure of the series. The Borel transform of 
Ro{Q^), satisfying 

RoiQ')= / 5(z)e-^/'^(«')rfz, (3) 

Jo 

as that of the observable R{Q'^), has singularities ("renormalons") on the real axis at 

z = p/Po, where p are integers, which are responsible for the factorial increase of the 

standard perturbative coefficients. 

On the other hand, a conformal expansion is totally free of renormalons. It 
may well contain other types of factorially increasing coefficients, e.g. due to the 
multiplicity of diagrams, but it cannot be influenced by running-coupling effects since 
the coupling in the conformal theory does not run. 

A simple demonstration of how conformal relations become free of renormalons 
is the following: consider (|l]) in the case that the coupling does not run, a(fc^) = Opp: 



the integral for i?o in (HD simply reduces to Opp (where the normalization of is 1 
by definition). If instead the coupling a(A;^) has an infrared fixed-point, changing 
variables in (0), e = fc^/Q^, one obtains 

Ro{Q')=l a(eQ^)0(6)-; (4) 

taking the limit Q^ — > inside the integral one obtains again R^^ = Opp, i.e. a 
trivial expansion. We shall return to this example below. 

The fact that conformal relations are renormalon-free becomes relevant to QCD if 
an all-orders relation exists between the standard QCD perturbative expansion and a 
conformal expansion. Indeed, such a relation exists when a non-trivial infrared fixed- 
point is realized perturbatively, as occurs [|ll]-|]l5l for sufficiently small Pq. Using the 



Banks-Zaks expansion, an A^j-independent conformal relation between two generic 
effective charges can be written [|ll| (see also [|ll]^[ll9l). The coefficients in such 



a relation are renormalon-free. On the other hand, these coefficients are explicitly 
expressed as combinations of the standard (non-conformal) perturbative coefficients 
with those of the (3 function. It immediately follows that in these combinations 
renormalon factorials must conspire to cancel out. In this way, the absence of renor- 
malons in conformal expansions translates into a constraint on the structure of the 
non-conformal perturbative expansion. 

After briefly recalling the relation between conformal expansions and the standard 
perturbative expansion, we demonstrate that a generic model for the Borel function 
satisfying the expected large-order behaviour of the perturbative series is not nec- 
essarily consistent with the conformal limit. Next, since the conformal constraint 
holds by definition in (|l|) we analyse this example in detail, showing explicitly how 
the renormalon factorials eventually conspire to cancel out, leaving the conformal 
expansion renormalon-free. 

Suppose that the perturbative expansion of R{Q^) is given by 

RiQ^) = aiQ^) + naiQ^f + r^aiQ'f + • ■ • , (5) 

where a{Q'^) satisfies a renormalization group equation (0) which has a non-trivial 
infrared fixed-point Opp = a{Q'^ = 0) = — /3o//3i + 0(/3o), for sufficiently small (3q (and 
/?i < 0). One can then write a conformal relation between the fixed-point values of 
R{Q^ = 0) = i?pp and app, 

i?pp = Opp + ciOpp^ + C2app^ + csOpp"^ + ■ • • , (6) 

The conformal coefficients q can be related to Tj in (|^) through the Banks-Zaks 
expansion 

Opp = ao + Vial + V2al + vso^ + ■ ■ ■ (7) 



where ao = ^-g-p — , and Vi depend on the coefficients of (3{a) [|1^, |I^ |18|. Since r, 
are polynomials of order i in Nf , one can write 

i 

'r'i = J2'^i,i^^o'', (8) 

k=0 



yielding the following A^j-independent relations fTT 



Cl = 


= ^^1,0 


C2 = 


= ri^i + r2,o 


C3 = 


= -ri^iVi + ^2,1 + ^"3,0 


C4 = 


= 2ri,iWi^ - ri^i^a - r2,if i + r2,2 + rg^i + r4,o. 



(9) 



For simplicity we shall assume in the following that a(Q^) satisfies the two-loop 
renormalization group equation. In addition, we shall ignore the Nf dependence of (3i. 
Under these assumptions Opp = Oq = — /?o/A) i-6- Wj = for any i >1. It obviously 
follows that 

[i/2] 
fc=0 

for any i, where the square brackets indicate a (truncated) integer value. In this model 
then, Ci is simply the sum of all the possible Vj^k coefficients such that j + k = i and 
j > k. Relaxing these restrictions is possible, but the price will be more complicated 
expressions. The conclusions would not change, provided that /3(a) itself does not 
contain renormalons (this is probably true in MS). 

Relation ([To| ) is intriguing: the q on the l.h.s. are conformal coefficients, which 
must be totally free of renormalons, while ri-k,k on the r.h.s are factorially increasing 
because of renormalons. The only way in which the condition that q are renormalon- 
free can be realized is if explicit cancellations occur on the r.h.s. To show that 
this cancellation is non-trivial, we shall now consider simple models for the Borel 
transform B{z), which are consistent with the expected large-order behaviour of the 
series and demonstrate that the corresponding q in ( p!OD become factorially increasing, 
thus violating the conformal constraint. 

We begin with the simplest example corresponding to a single simple pole in the 
Borel transform of the observable R{Q^): 

where Zp = p/ (3q is the renormalon location. Note that in this example we choose the 
renormalon residue to be a constant, but in fact in QCD it depends on Nf. 
The inverse Borel transform (§) yields 

R{Q^) = -ZpEiil, -^p/a)e-^-/^ (12) 

where a = a{Q^). The perturbative series of R{Q^) has the following factorially 
increasing coefficients 

(13) 

Under the assumption that (5i is A^/-independent, the decomposition of the coef- 
ficients of (|13D in powers of oq according to (H) yields Tj j = (— /5i/p)*i! and r, .,• = 




for any j ^ i. The resulting coefficients (0) in this model are therefore 

""''X («/2)! (-A/p)'/' I even ' ^ ^ 

Thus, the would-be "conformal coefficients" do diverge factorially. In some sense the 
factorial divergence is slowed down: q contains just (V2)! rather than %\. Conse- 
quently we define u = i/2 and write the Banks-Zaks expansion as: 

i?pp = ao E«! H^ ao'" = ao ^^K-^)^", (15) 

u=0 \ P / u=0 



where 

S = -w- (16) 

Po 



We found that in the simple Borel pole example the factorial divergence of the per- 
turbative series does enter the "conformal coefficients" . The conformal constraint is 
therefore explicitly violated. 

However, this example is not completely self-consistent: on one hand it was as- 
sumed that a((5^) runs according to the two-loop P function (it has a fixed-point), 
but on the other hand we used the large-A^/ (i.e. one-loop (3 function) form of the 
Borel singularity, namely a simple pole. It is known that a non-vanishing two-loop 
coefficient in the f3 function modifies the Borel singularity to be a branch point. For 
instance, for the leading infrared renormalon associated with the gluon condensate 
[p = 2) one has the following singularity structure in the Borel plane [Q 

where 6 is defined in (0). The corresponding perturbative coefficients are 



_ T{l + 5 + t) //?o' 

'' r(i + <5) U. 

The large-order behaviour is 



(1* 



which is different from ([TB|). 

As opposed to the previous example, the rj are not polynomials in Pq, so starting 
with ( P^ ) we cannot obtain the form (^. To see this, let us examine the expansion 
of the r function in r^ 

M6)^^,±^ = {6 + z){6 + z-l){6 + z-2)...{6 + l). (20) 



It is clear that fi{S) can be written as a sum 



MS) 



E /f ^'. 



(21) 



fc=0 



where fl are numbers. Since 5 is proportional to l//9o^, fiiS) contains all the even 



powers of l//?o from up to 2i. The additional positive power of Pq in (|Tq) finally 
leads to having half of the terms with positive power of f3o and half with negative 
powers. The latter correspond to non-polynomial functions of Nf, which cannot be 
obtained in a perturbative calculation. This suggests that the current example is 
unrealistic. 

The actual QCD situation, where the coefficients are polynomials in (3q and be- 
have asymptotically as ([T9|), can be imitated by truncating the negative powers of Po 
in ([T8|) . We verified explicitly (see below) that this truncation does not alter the 
eventual large-order behaviour of rj. Note that there is some ambiguity in the trun- 



cation: one can, in principle, truncate (pT[) at different k values and still obtain the 
same asymptotic behaviour. We shall choose the most natural possibility: truncate 
just the terms that lead to negative powers of Pq. 

In order to proceed we should find the coefficients fl . This can be done by 
writing a recursion relation using the property fi{6) = {6 + i)fi-i{6). This condition 
is equivalent to the following 



/f 



fk 



1 



^n 



(t-i) 



k = i 

< k <i 

k = 



(22) 



It is straightforward to use these recursion relations to obtain /^ to arbitrarily high 
order. 

After truncating the terms with negative powers of /5o, the coefficients (0) become 



r,; 




E /f s'- 

k=0 



(23) 



In order to verify that the truncation of the high powers of 6 does not affect the 
large-order behaviour of the series, we calculated the ratio 



ri/ri 



[[V2] _ 1 

E ft^s' 


1 


Y.ft's' 


_k=0 


La;=0 J 



(24) 



for various values of 5, as a function of the order i. It turns out that this ratio 
approaches 1 fast, indicating a common asymptotic behaviour. For instance, for 



b = 462/625, corresponding to eq. (|16D with p = 1 and the values of Po and Pi in 
QCD with Nf = 4, we find fi/ri ^ 0.995 for i = 8. 



Next we write the decomposition of fj as a polynomial in oq according to 



"^ij Jtii\ ^) 






(25) 



where j is odd for odd i and even for even i (as always, j < i). Finally, the "conformal 
coefficients" ( [T0| ) in this example are 



C2« 



j=0 



2m 



7 . fk-ui 1 



and the expansion is 



-Rf 



n=0 



2n 



.k=u 



E /ru-i 



\u+fc 



■k=u 






;-5)-". 



(26) 



(27) 



The square brackets should be compared with u\ in eq. (pISj), characterizing the simple 
Borel pole example. It turns out that the C2u increase /aster than u\, but slower than 
(2m)!. Thus the factorial behaviour of the "conformal coefficients" persists also in 
this example. Again, the conformal constraint is violated. 

Another possible approach to the analysis of the Borel cut example ([TTp is the 
following: consider the large-order behaviour of the coefficients in eq. (plQ]). Let us 
now ignore the l/r(l+5) factor, which does not depend on i and can be absorbed into 
the residue of the renormalon and expand i^ ~ exp(5 ln(i)) = 1 + 5 ln(z) + ^5"^ In^(z) + 
■ ■ ■. Again we find that high powers of 6 lead to non-polynomial dependence of the 
coefficients. As before we truncate these terms and write an approximation to rf*: 



^\tl}A^)S'. 



P 



{21 



k=0 



We checked numerically that the ratio fi/rf'^ approaches 1 as i increases, so the 
asymptotic behaviour is not altered by this truncation. Note that the powers of ln{i) 
in (p8D can be understood diagrammatically, as explained in 0. 
We proceed and write 



«j 



2 ■ 



-lyiimf-^ (^ 



i+j 
2 



p 



(29) 



and finally, using (0), the "conformal coefficients" are 

k\ 



C2u 



0=0 



■J,3 



2u 



E 



uik-u) 



-l)'=(lnA;)* 



'Pi 
.P . 



(30) 



The large-order behaviour of C2u turns out to be again between u\ and 2m!. In fact, 
the two ways we used to construct the coefficients in this example lead to roughly 



6 



the same asymptotic behaviour of the "conformal coefficients": the ratio between 
C2u and C2u approaches a geometrical series at large orders. The reason for this is 
simply the fact that in both examples the dominant term in the sum is the one 
coming from the highest power of the coupling (f2n,o in ( p6D and f2u,o in (pO|)). In 
fact, the contributions to the "conformal coefficients" from increasing orders in the 
coupling are monotonically increasing in both cases. We stress, however, that the 
decomposition of r^ into polynomials in /3o (^ij) is not at all similar in the two cases. 

We saw that, in general, consistency with the large-order behaviour of the pertur- 
bative series does not guarantee consistency with the conformal limit. In the above 
examples the renormalon factorials in Vij do not cancel out in the sum (|l^) but rather 
induce a non-physical factorial increase in the "conformal coefficients" . Apparently, 
the conformal constraint is non-trivial. 

The examples above also teach us that the large-order behaviour of the perturba- 
tive coefficients Tj (or the nature of the Borel singularity) by itself does not uniquely 
determine the decomposition of Tj into powers of (3q: several choices of Tj .,• can ffi. 
Indeed, we shall see below that the cancellation of renormalons in (0) crucially de- 
pends on having an "appropriate" decomposition of Tj. Since the model (|l|) is, by 
definition, consistent with the conformal constraint, the factorials should cancel out 



in (|10D and the model should give an example of an "appropriate" decomposition. 

In order to analyse the conformal coefficients in ([1|) we restrict ourselves to the 
contribution to Rq{Q'^) from small /c^, which is the origin of infrared renormalons, 
and expand the momentum distribution function 



{kVQ' 




(31) 



where $„ are numbers. It is sufficient to consider a specific infrared renormalon with 
n = p, so we choose our "observable" as 

where the upper integration limit is set for simplicity to Q^ . 

It was shown in pO|, ^ that if a(fc^) satisfies the two-loop renormalization group 



equation, the Borel representation of i?o is 

MQ') = l e-^'——--^e-^/^dz, (33) 

■^0 [1 - {z/Zp)] 

where a = a(Q^). This integral resums all those terms in eq. (|l]) that depend only on 
the ffist two coefficients of the (3 function. Note that eq. ( ^2]) is well defined thanks to 
the infrared fixed-point of the coupling a(fc^). On the other hand, eq. (|33| ) is not well 
defined because of the infrared renormalon, and it differs [^, ^ from eq. (^) by an 
ambiguous power correction. The equality between (|32|) and ( ^3]) should therefore be 
understood just as an equahty of the (all-order) power series expansion of the two 
expressions. 

7 



To expand (p3D we note that the Borel transform of Rq with respect to the 
modified coupling a, 

i^i + |. (34) 

coincides with the Borel transform (0). Using the coefficients (|18D we have 

"'-h' "S r(i + ^) Uj" • <'" 

Substituting a*"^^ for 



we obtain 






T(l + S + i) (!+*:)! /'/J„\7 A, ,^.t 



Defining n = k + i and performing first the summation over i we obtain 

oo 



n=0 



with the perturbative coefficients r„ given by 

^■.-(^)"i: ., '"./'a!:^ -^r-. (39) 

\Pj f^o^^■{n-^y■ r(l + 5) 

Note that the model eq. (|1^) corresponds to keeping only the i = n term in eq. (|39|). 
We now use (pT|) to expand the F function and write explicitly the dependence on Pq 
(or Oo). Defining j = 2i — 2k — n, we obtain 

n 

'^n= J2 '^njao^ (40) 

j=-n 

with 



-. = f§Vt7-4^(-i)7«.. (41) 



r, 



where u = ^^^. Next, we note that for < 2m < n, 

n! 



l=U 



n — i)\i\ 



E^^VT(-Wr. = 0, (42) 



so the negative powers of Oq are absent in eq. (|y). We thus identify a major difference 
between this example and the examples considered above: here the decomposition of 



Tn into powers of qq does not lead to any non-polynomial dependence, and truncation 
is not required. Note, however, that there are non-trivial cancellations. 

Finally, using (^), the conformal coefficients^ corresponding to Rq are given by 



2u 



C2n — 2^ 1^2u-j,j — 2^ rk. 



2u-k 



2u k 



k\ 



7 y 7 y // -Ni •] \ / Ji—u 



(i) 




(43) 



/ ^ ' ZU—JJ / , ' K,ZU — K / , / , /r 'M 'l 

i=0 k=u U=uJ=u l'^ ~ "*/•*• 

where the last equality was checked explicitly. In other words, the final result is 

i?o(Q' = 0) = app, (44) 

in accordance with our expectations. As explained above, the vanishing of the confor- 
mal coefficients in this case can be understood directly from the defining integral i?o 
(the polynomial N f dependence of the r„'s is also transparent from this representa- 
tion). We note that contrary to the previous examples (^) and (|30D , in (|43|) the 
term originating from the highest power of the coupling does not dominate. This is 
crucial for the eventual cancellation. 

In conclusion, we saw that the absence of renormalons in conformal coefficients 
can be seen as a constraint on the form of the perturbative expansion: renormalon 
factorials must cancel out in certain combinations. We considered various examples 
for the Borel transform, which are consistent with the same large-order behaviour, 
finding that this cancellation is non-trivial. It would be interesting to find a concrete 
general form of this constraint, which still appears elusive. We note that in the two- 
loop example studied here the constraint is implemented through the regular factor 
exp \ — Q^) in eq. (p^), and therefore it involves an infinite series of sub-leading terms 
in the large-order asymptotic behaviour of the coefficients. 
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